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Abstract 

The standard ACDM paradigm is complemented with a magnetized contribution whose 
effects on the anisotropies of the Cosmic Microwave Background (CMB) are assessed by 
means of a dedicated numerical approach. The accuracy on the temperature and polarization 
correlations stems from the inclusion of the large-scale magnetic fields both at the level of the 
initial conditions and at the level of the Einstein-Boltzmann hierarchy which is consistently 
embedded in a generalized magnetohydrodynamical framework. Examples of the calculations 
of the temperature and polarization angular power spectra are illustrated and discussed. The 
reported results and the described numerical tools set the ground for a consistent inclusion 
of a magnetized contribution in current strategies of cosmological parameter estimation. 



Current analyses of cosmological data sets (see, for instance, [1]) are customarily per- 
formed using a variety of theoretical models which represent diverse completions (i.e. delicate 
improvements) of the pivotal A cold dark matter paradigm (ACDM in what follows). As an 
example, if we ought to know how large could be the contribution of a stochastic background 
of gravitational waves to the anisotropies of the Cosmic Microwave Background (CMB in 
what follows), the possible presence of tensor modes should be added, as a supplementary 
feature, to the basic list of cosmological parameters of the standard ACDM paradigm whose 
updated version (improved by a tensor contribution) will then be compared with the ex- 
perimental data. From the latter comparison likely values of the ratio between the tensor 
and scalar power spectra can be inferred. Other possible completions of the ACDM lore 
might include, for instance, massive neutrinos, a minimal duration of the inflationary phase, 
an effective barotropic index for the dark energy component and many others. Yet another 
class of delicate improvements of the ACDM paradigm contemplates the inclusion of one (or 
more) non-adiabatic modes which can be either correlated or anticorrelated with the stan- 
dard adiabatic component (see, for instance, [2]). For specific choices of the non-adiabatic 
amplitude and spectral index the fit to the experimental data may even improve (see, for 
instance, last reference in [2]). 

In this paper we are going to describe another completion of the ACDM paradigm. The 
possibility we are going to present will be dubbed as magnetized ACDM scenario (mACDM 
in what follows). Large-scale magnetic fields arise over different scales ranging from galaxies 
to clusters [3]. Superclusters have also been claimed to have magnetic fields [3] at the piG 
level even if firmer evidence is still lacking. Hopefully some of the findings of the Auger 
project [4] could be used for a magnetic "tomography" of the local Universe, say within a 
cocoon of 60 Mpc. 

If the present magnetized structures emanate from pre-existing cosmological relics they 
must be present prior to matter-radiation equality, affecting, in this way, the physics of 
photon decoupling and, ultimately, the formation of the CMB anisotropies. In the mACDM 
scenario the CMB anisotropies are computed in the presence of a stochastic magnetic field 
1 which will affect both the initial conditions and the dynamical evolution of the Einstein- 
Boltzmann hierarchy: 

(Bi&Bjip)) = ^^(k+^P^P^k), 

^■(*) = <Sy-^ PB(k) = A B (±) . (1) 

The minimal mACDM scenario contains, on top of the (six) parameters of the ACDM 
paradigm, two new parameters, namely the magnetic spectral index ub and the ampli- 

1 A stochastic magnetic field does not break, by definition, the overall spatial isotropy of the background 
geometry and allows, as a consequence, for angular power spectra just expressed in terms of their dependence 
upon the multipole I without any further preferred direction. For the opposite situation see, for instance, 
[5]. 
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tude of the magnetic power spectrum Pb(&) evaluated at the magnetic pivot scale k^. Non 
minimal extensions of the mACDM scenario include: the correlation (or the anticorrelation) 
between the adiabatic mode and the magnetic field, the simultaneous presence of a mag- 
netized adiabatic mode together with one (or more) magnetized isocurvature modes. For 
sake of simplicity we will stick here to the minimal situation reporting elsewhere on the 
complementary cases. There have been lately theoretical and semi-analytical works along 
this direction [6] but the moment has now come to tailor, for the first time, a consistent 
numerical approach to treat the effects of a fully inhomogeneous magnetic field on the scalar 
modes of the geometry. 

Since, after neutrino decoupling, the concentration of electrons and protons is roughly 10 
orders of magnitude smaller than the concentration of the photons, the Debye length-scale 
will be approximately 20 orders of magnitude smaller than the Hubble radius (for instance at 
matter-radiation equality). The conductivity of the globally neutral plasma, i.e. <j c , will be 
dominated by Coulomb scattering so that a em a c ~ 47r(T/m c ) 1 / 2 T(ln A c ) _1 where A c is the 
argument of the Coulomb logarithm. Ohmic electric fields are then suppressed with respect 
to the total Ohmic current, as it is the case in a good terrestrial conductor: 

* - ft J V x B 

E + v h x B = — ~ , a = aa c , (2) 

a Ana 

— » _ — * — * ( , — » 

where E = a £ and B = a 2 B are the electric and magnetic fields rescaled through the second 
power of the scale factor a(r) of a (conformally flat) Friedmann- Robert son- Walker (FRW) 
geometry. The large-scale description of our problem is, physically, the curved-space version 
of MHD where the electric field, the magnetic field and the Ohmic current are all solenoidal 
quantities. The effects of the magnetic field on the scalar modes of the geometry are far 
from trivial: the evolution of the bulk velocity of the plasma is given by 2 

0b + 70b = ^V(fty ~ ft,) + V ' [ 4 e h = V-V h . (3) 

6 pb a Ph 

Given the hierarchy between the electron and proton masses the bulk velocity of the plasma 
#b will be essentially the proton velocity. Because of the Coulomb-dominated conductivity, 
the magnetic fields will be present and not diffused for typical length-scales larger than the 
magnetic diffusivity length L a whose ratio to the Hubble radius is so small (i.e. HL a ~ 
3.9 x 10 _17 (T/eV) 1 / 4 ) that the k a ~ L^ 1 will provide effectively an ultra-violet cut-off to 
the magnetic power spectrum of Eq. (1). In the presence of diffusion damping (i.e. shear 
viscosity) the Silk wave-number will be the dominant dissipative scale for the baryon-photon 
fluid. Under these conditions the magnetic flux (and helicity) will be effectively conserved 
to a very good approximation. 

2 A prime denotes a derivation with respect to the conformal time coordinate r and TL = (In a)'. The 
notation 6x = V • vx denotes the three-divergence of the peculiar velocity of the species X. Throughout 
the paper e' = x c n c <jThn c a/ao denotes the optical depth. Finally, with similar notations, 5x = S s px / px is 
the density contrast of the species X. 
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Large values of the conductivity break explicitly Lorentz invariance. The plasma frame, 
where the electric fields are suppressed in comparison with the magnetic fields, arises natu- 
rally. The scalar fluctuations of the geometry, still relativistic, will be treated, numerically, 
in the synchronous gauge (S'-gauge in what follows) where the only non-vanishing entries of 
the perturbed metric are, in Fourier space: 



S s gij(k, r) = a 2 (r) kikjh(k, r) + 6£(fc, r) (kikj - -Sij 



(4) 



where the symbol S s emphasizes that we are dealing here with scalar (as opposed to vector or 
tensor [7]) fluctuations of the geometry. The mACDM code is an extension of the CMBFAST 
package [8] which is, in turn, based on the COSMICS package [9]. The evolution of the scale 
factor is integrated numerically from the usual Friedmann-Lemaitre equations 

^ = 8^ fl2pt) H 2 -H' = ±nGa 2 ( Pt + p t ), (5) 

where p t and p t are, respectively, the total energy density and pressure of the plasma. 
Equation (3) will then be supplemented by the governing equations for the magnetic fields 
as well as by the CDM particles and by the neutrino component: 

5' C = -9 C + ^, 6> C + W C = 0. (6) 

K = -\e v + 2 -h\ ff v = - iv 2 5„ (7) 

4 2 4 

a' = —ft, - —ti - (8) 
" 15 15 5 s " K ' 

where a u = J-v2/2 is the quadrupole of the perturbed phase space distribution. Higher mul- 
tipoles of the Boltzmann hierarchy (like the octupole J- V z) follow by setting initial conditions 
on the lower multipoles. The governing equations for baryons and photons are given by Eq. 
(3) together with 

0, = + e'(6 h - ft,), 5; = -U, + h>, 5' h = -6 h + ^+ E '£** ) , (9) 

where the last term in the evolution of 5^ is negligible at finite conductivity. The plasma 
and the magnetic fields all gravitate and contribute to the Hamiltonian and momentum 
constraints whose specific form is, respectively: 

2V 2 £ + Hh' = -8irGa 2 [5 sPt + 5p B ], (10) 

^ = ^Ga 2 {( Pt + Pt )9 t + ^^}, (11) 

where 5 s p t is the total density fluctuation in the S gauge and (p t + pt)&t — DaG°a + Pa)^a 
is the total peculiar velocity. Since the conductivity a is always large, the contribution of 
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the MHD Poynting vector is, in practice, always negligible. The (ij) components of the 
perturbed Einstein equations read 

h" + 2Hti + 2V 2 £ = 247rGa 2 [5 Pt + 5p B ], (12) 
(h + 60" + 2H(h + 60' + 2V 2 £ = 2ArrGa 2 [( Pu + p u )a u + (p 7 + p 7 )<7 B ], (13) 

where 5 s p t is the fluctuation of the total pressure. In MHD, AttJ = V x B so that the Lorentz 
force and the magnetic anisotropic stress (associated with a B ) are related by: 



-V • [(V x B)xB} + 



V 2 ft T 



5p B (r,x) 



(14) 



The simplest set of initial conditions to be imposed on the hierarchies of the fluctuations in 
the intensity and of the polarization is the magnetized adiabatic mode. When Coulomb and 
Thompson couplings are both tight (i.e. 9 h ~ # 7 = # 7b ) the whole system of the governing 
equations can be solved in the limit when the relevant wavelengths are larger than the Hubble 
radius prior to matter-radiation equality (i.e., in terms of the wave-number fc, ^t<1). To 
lowest order in kr the magnetized adiabatic mode reads, in Fourier space: 



h(k, t 



e c (k,r 

a u (k, t 
where R v = r / 
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r + 1) with r = 0.681 (iV„/3) (and i? 7 



R v ). The constant mode of h 
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Figure 1: The magnetized temperature autocorrelations stemming from the magnetized 
adiabatic mode (left plot) and in the absence of an adiabatic component (plot at the right). 
The other parameters are fixed to the central values of the best fit of WMAP three year 
data alone (plot at the left) and to the central values of the best fit of the WMAP data 
supplemented by the gold sample of supernovae (see [1]). 



leads to a gauge mode and should be projected out of the solution. The second potentially 
dangerous gauge mode is fixed by setting to the CDM peculiar velocity 6 C . The spectrum of 
C(k) is simply related to the spectrum of 7Z which parametrizes the curvature perturbations 
on comoving orthogonal hypersurfaces [6]. Deep in the radiation epoch, according to Eq. 
(15), TZ(k,r) = —2C(k). Consequently the initial conditions for the adiabatic component 
will be given in terms of the spectrum of TZ, i.e. V-ji(k) = Anik/k^) 11 ^ 1 where n s is the 
adiabatic spectral index and An is the amplitude of the scalar power spectrum at the pivot 
scale k p = 0.002 Mpc -1 . The spectrum of the magnetized contribution will be encoded in 
the spectra of <tb and i-e. respectively, V a {k) = 5(?iB)fi Bi (fc/fcL) 2 ' n °" 1 ' and Vn(k) = 
J r (nB)^BL(k/kh) 2 ^ nB l ^ where k^ denotes the magnetic pivot scale 3 which will be taken 
1 Mpc -1 . In the minimal mACDM scenario n a = n-Q. The spectra Vn(k) and V a (k) are not 
arbitrary but rather computed from Eq. (1) through a rather standard procedure (see, for 
details, [10]): 

n BL = ^ = 7.5xio- 9 (% 

87rp 7 VnG 

3 The range of magnetic spectral indices discussed here is 1 < «b < 5/2. In this case, regularizing the 
magnetic field with a Gaussian window function, the spectra of the 17b and <tb are the ones reported in Eq. 
(25). For values outside the mentioned range cut-offs are required either in the infra-red or in the ultra-violet 
[10]. The value of the magnetic pivot scale implies that £?l coincides effectively with the proper amplitude 
of the magnetic field regularized over a Mpc window at the onset of protogalactic collapse [6]. 
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Figure 2: The magnetized TE and EE correlations are illustrated for different values of the 
spectral indices and different values of the magnetic field intensities. 



„, N (2vr) 2 ("B-i) 
G(n B ) - 



r 2 
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n B + 29 



15(5-2n B )(n B - 1) 



20(7 -n B ) 
F(n B ) = nB + 29 Q{n B ). (25) 



In Fig. l(plot at the left) the effects of the magnetized contribution on the temperature 
autocorrelations is illustrated for different values of the spectral index and of the magnetic 
field intensity. If the strength of the magnetic field is augmented (or the value of the spectral 
index becomes bluer) a distortion of the second and of the third peaks is correlated with an 
increase of the first peak. This spectral distortion fits with the results of a semi-analytical 
calculation conducted in the longitudinal gauge (see [6], last reference). In spite of the correct 
spectral distortion on the shape, the semi-analytical argument is intrinsically less accurate. 
The predicted height of the first peak varies in a non-monotonic way with the variation of 
the spectral index and of the magnetic field intensity [10]. If the adiabatic mode is absent 
from the initial conditions (Fig. 1 plot at the right) the amplitude of the TT correlations 
is always much smaller than in the case of the magnetized adiabatic mode. Furthermore 
a hump appears at intermediate multipoles. In Fig. 2 we illustrate the angular power 
spectrum of the temperature-polarization cross-correlations (for short the TE correlations) 
as well as the angular power spectrum for the polarization autocorrelations (for short the 
EE correlations). According to Fig. 3, for high i (where, hopefully, there will be, in the near 
future more precise data from the Planck explorer mission) the TT correlation is definitely 
sensitive, for i ^> 1500, to a nG magnetic field. The TE correlation is also sensitive to the 
nG range (see Fig. 3, plot at the right) especially as soon as the spectral index increases 
from the nearly scale- invariant limit. The magnetized CMB anisotropies have never been 
computed consistently and systematically for the scalar modes of the geometry which are the 
ones observationally more relevant. In this paper we built a consistent numerical approach to 
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Figure 3: The TT and TE correlation for t ^> 1500. The full curve is the extrapolation of 
the three year best fit of the WMAP data alone. 

the problem and presented a magnetized completion of the ACDM paradigm. Our approach 
is accurate enough to resolve nG magnetic fields. The usual custom of setting a bound on 
the magnetic field becomes too simplistic in the light of the present results and according 
to the standards of CMB physics. In a similar perspective, for instance, the bound on the 
ratio between tensor and scalar power spectra can only stem from an appropriate strategy 
of parameter extraction where the tensor contribution is included in the CMB anisotropy 
calculations and fits. The same must be done in the case of the mACDM scenario. If the 
correct strategy is not enforced, potentially interesting degeneracies between the parameters 
of the magnetized background and other cosmological parameters can be totally overlooked. 
Consequently, our program is, in the short run, to extend and complete our code to the 
case of tensor and vector modes (which are known to be far less relevant at large scales) 
and to the case of the various magnetized isocurvature modes. In parallel we ought to test 
the mACDM scenario against the various data sets eagerly waiting for Planck data and its 
claimed high accuracy for large multipoles. 

K.E.K. is supported by the "Ramon y Cajal" program and grants FPA2005-04823 and 
FIS2006-05319 of the Spanish Science Ministry. 
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